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ABSTRACT 


Although sporadic meteoroids generally pose a much 
greater hazard to spacecraft than shower meteoroids, me- 
teor showers can significantly increase the risk of dam- 
age over short time periods. Because showers are brief, it 
is sometimes possible to mitigate the risk operationally, 
which requires accurate predictions of shower activity. 
NASA’s Meteoroid Environment Office (MEO) generates 
an annual meteor shower forecast that describes the vari- 
ations in the near-Earth meteoroid flux produced by me- 
teor showers, and presents the shower flux both in abso- 
lute terms and relative to the sporadic flux. The shower 
forecast incorporates model predictions of annual varia- 
tions in shower activity and quotes fluxes to several limit- 
ing particle kinetic energies. In this work, we describe 
our forecasting methods and present recent improve- 
ments to the temporal profiles based on flux measure- 
ments from the Canadian Meteor Orbit Radar (CMOR). 
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1. INTRODUCTION 


Meteoroid impacts are known to cause damage to space- 
craft surfaces; their high speed (12-72 km s~!) compared 
to orbital debris means that relatively small particles can 
be hazardous. Most of the meteoroid flux is associated 
with “sporadic” meteoroids, which are those not associ- 
ated with any meteor shower. The sporadic complex is 
present throughout the year and constitutes the vast ma- 
jority of microgram-or-larger meteoroids. For this rea- 
son, meteoroid environment models such as NASA’s Me- 
teoroid Engineering Model (MEM) [1, 2] and ESA’s In- 
terplanetary Meteoroid Environment Model (IMEM) [3] 
focus on the sporadic complex. Although these models 
include showers in the total meteoroid flux, they do not 
model the short-duration fluctuations caused by showers. 
Such fluctuations are generally not worth modeling for a 
long-duration spacecraft mission. 


Over short time spans, however, meteor showers can 
match or even exceed the sporadic flux. The Geminid 


meteor shower, for example, can double the meteoroid 
flux at the time of peak activity. The Leonid meteor 
shower does not pose much of an impact risk on a typical 
year, but has occasionally produced outbursts in which 
the number of meteors is tens to thousands of times larger 
than normal [4]. In such cases, operational spacecraft 
sometimes choose to mitigate the risk by avoiding op- 
erations that increase the spacecraft’s vulnerability. If the 
risk is great enough, operators may consider re-orienting 
the spacecraft to present its least vulnerable side to the 
shower, phasing its orbit to use the Earth as a shield, or 
powering down components. 


All mitigation strategies have some associated cost such 
as delays, reduced functionality, or fuel. Accurate shower 
predictions are therefore needed for informed risk assess- 
ment. The Meteoroid Environment Office (MEO) pro- 
duces annual meteor shower forecasts that can be used 
to assess the increase in the particle flux due to showers. 
These forecasts are designed to be used in conjunction 
with an existing sporadic meteoroid risk assessment; un- 
like sporadic models, however, the annual forecast pro- 
vides a detailed temporal profile that takes showers into 
account. The forecast also takes shower variability into 
account; any unusual activity predicted by modelers [5, 6] 
is incorporated. 


While modelers can often predict the level of meteor 
shower activity, we frequently rely on past observations 
to characterize the activity profile. For many years, the 
MEO (and others; [7]) used a set of activity profiles 
that were determined from naked-eye meteor observa- 
tions [8]. Many of the activity profiles in this set are quite 
noisy due to low number statistics. Furthermore, daytime 
meteor showers were of course not visible to observers 
and thus the MEO had to employ a default activity pro- 
file for many of these showers. 


Meanwhile, the Canadian Meteor Orbit Radar (CMOR) 
began operation in 2001 [9]. CMOR is a patrol radar 
and can detect meteor echoes during the day as well as at 
night. It is located near London, Ontario and thus moni- 
tors the northern hemisphere, although it can detect me- 
teors with declinations as low as -40°. Meteor flux mea- 
surements from CMOR [10] are available for the years 
2002-2015, providing us with a data set suitable for re- 
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Figure 1. Sample meteor shower activity profile. The 
peak zenithal hourly rate is ZHRo and the peak occurs at 
solar longitude \o. On either side the activity follows an 
exponential function of solar longitude; the two parame- 
ters By, and By, are not necessarily equal. 


vising many shower activity profiles. 


This paper presents both our methodology for shower 
forecasting (Section 2) and the improvements we’ve 
made to our shower activity profiles based on CMOR flux 
measurements (Section 3). Finally, Section 4 discusses 
the forecast and the changes in it from a spacecraft risk 
perspective. 


2. FORECASTING METHODS 


This section describes our algorithm for forecasting me- 
teor shower fluxes. These fluxes are derived from a set 
of meteor shower parameters and are adjusted for an al- 
titude of 400 km. We also compute the appropriate spo- 
radic meteor flux at this altitude, taking gravitational fo- 
cusing and planetary shielding into account. Finally, we 
compute the flux enhancement produced by these show- 
ers over the baseline meteoroid flux. 


2.1. Shower fluxes 


Meteor showers do not have a well-defined duration. In- 
stead, peak activity occurs at a particular time each year 
and activity gradually increases up to the peak and de- 
creases after the peak (see Fig. 1). A double exponential 
function has been found to be a good fit to most meteor 
shower activity profiles [8]. Measures of effective shower 
duration (such as the full width at half maximum or time 
period over which activity exceeds a given threshold) thus 
depend on the steepness of this double exponential func- 
tion. 


We characterize the activity profile of each shower with 
four parameters: peak time (Ao), peak zenithal hourly rate 


(ZHRg), and two exponents that characterize the shape 
of the shower’s activity profile (B, and B,,). The time 
of peak shower activity occurs when the center of the 
meteoroid stream intersects the Earth’s orbit; we there- 
fore measure time in terms of solar longitude, Aj. The 
zenithal hourly rate (ZHR; the rate at which meteors oc- 
cur when the shower radiant is directly overhead) in- 
creases with time before the peak and decreases after the 
peak: 
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ZHR = ZHRy - { (1) 


In some cases, such as the Perseids, the activity profile 
of a single meteor shower is constructed from two sets of 
parameters that describe “peak” and “base” activity. 


ZHR describes the rate at which a meteor shower pro- 
duces visible meteors. ZHR can be converted to meteor 
flux by taking into account observer biases and shower 
characteristics. We use the methodology of [11] to cal- 
culate the flux of meteoroids that have a brightness of at 
least magnitude 6.5: 


_ ZHRo - Kayg + (13.17 — 16.5)(r — 1.3)9-748 
37200 km? 
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where the average perception factor, kaye, is of order 
unity. The population index, r, describes the brightness 
distribution of meteors within the shower and is also re- 
quired as an input for each shower. Because ZHR has 
nls of hr—!, this equation yields flux in units of km~? 
hr-*. 


This magnitude-limited flux can be converted to a mass- 
limited (milligram or larger) flux as follows [11]: 


-1 
fie = tex 7 p75 logio (29 km s /v100) (3) 


This equation makes use of Verniani’s relationship [12] 
between magnitude, mass, and velocity to calculate the 
meteoroid mass that produces a magnitude 6.5 meteor at 
the shower’s speed. The shower velocity is that at the top 
of the atmosphere. 


Finally, the flux can be scaled to any arbitrary limiting 
mass using the relation: 
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where s = 1 + 2.3 log, r is the shower mass index. 


In all annual meteor shower forecasts to date, an addi- 
tional factor of 2 was applied to Eq. 2. This factor was 
obtained by applying Eq. 2 to the Griin interplanetary flux 
[13] and comparing it with a sporadic ZHR estimate of 8. 
A factor of 2 was found to bring them in rough agree- 
ment. However, a separate mass-luminosity relationship 
was used to convert magnitude to mass in that calcula- 
tion. When Eq. 3 is used instead, the flux corresponding 
to a sporadic ZHR is a factor of 4 larger, obviating the 


Table 1. The four kinetic energies (KE;e¢) to which the 
MEO annual meteor shower forecast reports fluxes. The 
second column lists the particle mass which, at 20 km 
s—1, has the listed kinetic energy. The third column lists 
the particle diameter which, for a bulk density of 1000 kg 
m—®, has the listed mass. 


KE vet J ) Meret (g) ret (cm) 
6.7 3.35 x 10-5 0.04 
104.7. 5.24 1074 0.1 
2827. 1.41 x 107? 0.3 
104720. 5.24 x 107! 1.0 


need for a factor of 2. Additionally, the meteoroid speed 
distribution [1, 14, 15] is much broader than the Griin 
assumption of a constant 20 km s~!. When we fold a 
speed distribution through Eq. 2, we obtain a number that 
can differ from that obtained using a single speed of 20 
km s~? by a factor of 2. Finally, the ZHR of ~ 8 that 
was used for this calibration may originate from [8] and 
is thus based on an hourly rate (HR), not ZHR. A later 
sporadic ZHR estimate by Rendtel took sporadic source 
radiants into account and yielded ZHR~ 22 for the apex 
source alone.[16, 17] For these reasons, we conclude that 
the sporadic ZHR does not present a robust calibration 
point for Eq. 2 and therefore drop the factor of 2 that was 
previously applied. 


Most ballistic limit equations (BLEs) are proportional or 
near-proportional to kinetic energy, not mass [18]. There- 
fore, the MEO shower forecast quotes fluxes correspond- 
ing to four characteristic kinetic energies, listed in Table 
1. Put another way, shower fluxes are calculated for four 
limiting meteoroid masses which relate to our reference 
masses in the following way: 
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Mim = Mref X ( 


where v490 is the shower’s speed at an altitude of 400 km. 


Historically, the meteor shower forecast has applied an 
“obscuration factor” of 0.6 (corresponding to an altitude 
of 400 km) to each shower flux to account for the fact 
that the Earth can, at times, shield the spacecraft from the 
shower. We omit this term from future forecasts for two 
reasons. First, this shielding effect is velocity-dependent 
[19], and therefore a constant value of 0.6, while close 
to the correct values, is not strictly correct. Second, it is 
important to provide spacecraft operators with a “worst- 
case scenario;” such a scenario would involve full expo- 
sure to the meteoroid stream at the time of peak activity. 
(This obscuration factor did largely negate the factor of 2 
that was previously applied to Eq. 2, resulting in shower 
fluxes that were over-conservative by only 20%.) 


We do include gravitational focusing, wherein the flux 
of slow showers in particular is increased at the top of 
the atmosphere due to Earth’s gravitational pull. Because 
shower fluxes are measured at the top of the atmosphere, 


we must invert the gravitational focusing effect to obtain 
the slightly lower shower flux corresponding to an alti- 
tude of 400 km [19]: 
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This effect is small: we find that for our slowest shower 
(the Draconids, at 20 km s~!), the flux at 400 km altitude 
is 98.6% of the flux at 100 km altitude. This difference is 
within our uncertainties but we include it for the sake of 
completeness. 


In addition to computing fluxes, the meteor shower fore- 
cast code integrates each flux to calculate the meteoroid 
fluence over a given time interval. This time interval is ei- 
ther the default value of 7 hours or an interval requested 
by a customer. 


2.2. Sporadic flux 


In order to facilitate risk assessments, the meteor shower 
forecast compares the total shower flux at a given time 
with the total time-averaged meteoroid flux. The total 
time-averaged flux includes both shower and sporadic 
meteoroids; as a result, periods of low shower activity 
correspond to a reduction in risk, while periods of high 
shower activity correspond to an increase in risk. 


The total meteoroid flux in interplanetary space is cal- 
culated using Eq. A3 of [13], which was derived assum- 
ing an interplanetary meteoroid speed distribution with 
Ums = 20 km s~t. An interplanetary speed of 20 km 
s~' corresponds to a speed of 22.75 km s~! at a lo- 
cation 400 km above the Earth’s surface. We therefore 
calculate four reference sporadic masses using Eq. 5 and 
V400 km = 22.75 km gt, 


We convert the Griin interplanetary flux to one valid 
at 400 km altitude by taking gravitational focusing and 
planetary shielding into account. The enhancement of 
the interplanetary flux caused by gravitational focusing 


is given by [19]: 
2 


V400 
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where v499 = 22.75 km s~! is the speed at 400 km and 
vp = 20 km s~! is the interplanetary meteoroid speed. 
The fraction of the interplanetary flux that remains after 


planetary shielding is given by: 
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where v is meteoroid speed at the top of the atmosphere 
and Rg = 6371 km is the radius of the Earth. 


The resulting sporadic fluxes for our chosen altitude of 
400 km range from 4.18 x 107? m7? yr7! for our small- 
est particle energy to 1.93 x 1077 m7? yr7! for our 
largest particle energy. 


Previous versions of the shower forecast used Eqns. 7-5 
and 7-6 of [20] to compute the gravitational focusing and 
shielding at 400 km. However, Eq. 7-6 is not correct for 
a interplanetary meteoroid speed of 20 km s~!, nor is it 
correct for the speed distribution presented in [20]. Addi- 
tionally, the forecast previously did not use the 22.75 km 
s~! speed at 400 km to adjust the reference masses; once 
again, these two choices fortuitously had opposite effects, 
resulting in a baseline flux that was close to accurate. 


2.3. Enhancement factors 


Finally, the meteor shower forecast reports flux enhance- 
ment factors, which quantify the increase in flux pro- 
duced by showers over the baseline Griin flux. We as- 
sume that the Griin flux includes showers, but is time- 
averaged, while the shower flux varies with time. Addi- 
tionally, Griin’s flux equation [13] reports the meteoroid 
flux incident on a randomly tumbling flat plate in inter- 
planetary space. The shower fluxes discussed in Section 
2.1 represent those on a plate facing the shower radi- 
ant; the zenithal hourly rate applies when the radiant is 
directly above the collecting surface. This combination 
produces a worst-case scenario shower-to-sporadic flux 
ratio, in which the spacecraft surface faces the meteor 
shower and does not benefit from planetary shielding. 


Both the orientation requirement and lack of shielding 
must be discarded, however, in order to calculate the av- 
erage shower contribution to the total meteoroid flux. Let 
us assume that the Griin flux (fg) consists of a sporadic 
background flux (fp) plus a time- and direction-averaged 
shower flux (f;). The intercepted shower flux is propor- 
tional to cos 0, where 6 is the angle between a surface’s 
normal vector and the shower radiant. If we average cos 0 
over all directions, and compare to the flux intercepted 
when @ = 0, we obtain the following ratio: 
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If we also include the reduction in flux produced by plan- 
etary shielding (7,), the fraction of the Griin flux made up 
of shower meteors is a = f,/ fo = +(f400(Ao)ns)/ fo; if 
faoo represents the shower flux on an un-shielded plate 
facing the shower radiant at 400 km altitude. This in turn 
allows us to calculate the sporadic flux as fp = (1—a) fa. 


Finally, the fractional enhancement of the flux due to 
showers at any point in time is: 


— faool(Ao) + fo - fe — favo(Ao) — of 
g ; 
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The forecast reports this enhancement factor as a per- 
centage of the Griin flux. Note that we use f4o9, rather 
than 7 faoons, to calculate g. The enhancement factors 
therefore reflect the increase in meteoroid flux on an un- 
shielded plate facing the shower radiant, which is in keep- 
ing with our goal of providing worst-case shower fluxes. 


(1) 


We calculate a using a set of standard meteor shower 
parameters that represent shower activity in a “typical” 
year. However, in any given year, individual showers 
could exceed or fall short of the standard, in which case 


+(f400(Ao)ns)/ fo # a. 


The meteor shower forecast reports the total shower en- 
hancement as a function of time, taking all active showers 
into account. Our g is thus technically an overestimate, 
as a flat surface cannot be normal to multiple radiants at 
once. However, one shower usually dominates at any 
given time and g therefore remains a useful estimate of 
the increase in meteoroid flux. The meteor shower fore- 
cast also provides radiant information for the most sig- 
nificant showers. If the flux enhancement exceeds a mis- 
sion’s risk tolerance, the forecast fluxes can be used in 
combination with these radiants to refine the meteoroid 
impact risk assessment. 


The meteor shower forecast has historically assumed that 
the Gr un flux contains both sporadic and shower meteors 
[13]. This choice was made based on personal commu- 
nication with the authors of that work, who claimed to 
have made no attempt to separate the two types of mete- 
ors. We plan to revisit this assumption in the future and 
determine whether the Griin flux equation is consistent 
with or exceeds sporadic-only flux measurements. 


3. FLUX MEASUREMENTS AND ACTIVITY 
PROFILES 


The lists of shower inputs used to generate MEO meteor 
shower forecasts are derived from a variety of sources. 
The shape parameters (Ag, By, and B,,,) for most show- 
ers are taken from [8], who fit double-exponential activ- 
ity curves to visual observations of meteor showers. His 
data were provided by observers in both the northern and 
southern hemispheres. However, visual observers are not 
able to measure rates for daytime showers unless they oc- 
cur near dawn or dusk, and are in general less able to 
measure radiants accurately. 


In this work, we use 14 years of CMOR fluxes to revise 
our meteor shower activity profiles. These fluxes are sin- 
gle station, but the sporadic contamination has been sub- 
stantially reduced by placing a detection angle restriction 
on meteor detections. Radar detection of a meteor re- 
quires that the meteor trail be perpendicular to the line- 
of-sight. As a result, meteors associated with a particu- 
lar radiant are detected along a great circle that lies in a 
plane perpendicular to that radiant. Selecting only those 
meteors that lie along this great circle removes a large 
number of sporadic meteors but preserves shower mete- 
ors. When multiple showers are active, these circles can 
intersect; such points of intersection have been excluded. 
Despite these measures, some background contamination 
persists, which we remove in Section 3.1. 


CMOR’s gain pattern and radiant coverage have been 
taken into account to compute fluxes; see [10] for addi- 
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Figure 2. Raw fluxes corresponding to the AAN radiant 
measured by CMOR. The data are color-coded by year, 
where early years are blue and later years are red. Note 
the single, light-green outlier with a flux of nearly zero. 


tional details. Some uncertainty remains regarding the 
magnitude of these fluxes, but relative fluxes are thought 
to be robust and appropriate for constraining the shape, if 
not the peak value, of shower activity profiles. 


CMOR provides average fluxes for one-day intervals, and 
so a single year of data cannot provide a tight constraint 
on the exact peak time of the shower. However, the so- 
lar longitude corresponding to a given day varies by year, 
and so by combining many years of data, we can in some 
cases measure the peak time with less than 1-day resolu- 
tion. 


3.1. Detrending 


Fig. 2 shows the raw flux over all years for an example 
shower: the a Antliids (AAN). Several issues with the 
data are apparent: [1] the shower appears as an enhance- 
ment over the background flux, [2] the background flux 
appears to contain a trend with solar longitude, [3] the 
background flux tends to decrease over the years, and [4] 
anomalous outliers are present in the measurements. 


Our first task is to characterize and remove the back- 
ground flux, to which we fit a simple linear trend. In 
order to minimize the shower’s influence on this fit, we 
apply the following weights, w to the data: 


w= (Ao = ee (12) 


where Aj is solar longitude and Apom is the solar longi- 
tude corresponding to an estimate of the peak time. We 
then subtract this linear trend from the data. We find that 
this simple approach works fairly well to remove most of 
the background flux and nullify trends within it. 


We refine our removal of the background flux by fitting 
for an offset between the average background flux and 
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Figure 3. Detrended fluxes corresponding to the AAN ra- 
diant measured by CMOR. The data are color-coded by 
year, where early years are blue and later years are red. 


0.015 


the background flux per year. In this case, we attempt 
to exclude the shower and outliers by obtaining the aver- 
age flux for values between the 10th and 50th percentiles. 
Thus, if a shower spans more than half of the observation 
period, the background flux will be overestimated. Fig. 3 
shows the result of our background removal approach. 


3.2. Fitting 


We next fit a double-exponential profile (see Eq. 1) to 
our detrended data using y?-minimization. We fit for 
the peak time, peak flux, and the two exponents govern- 
ing the rise and decay of the shower (B+ and B-). Be- 
cause CMOR measures flux over day-long periods, we 
also convolve this double-exponential profile over a one- 
day period before fitting it to the data. A fit to another ex- 
ample shower, the Daytime Sextantids (DSX), is shown 
in Fig. 4. 


In our example as well as for a number of other showers, 
the shower amplitude appears to vary with year. This may 
be related to the decrease in the background over time. 
When this trend is apparent in a set of shower fluxes, we 
fit the amplitude of the double-exponential profile to each 
year’s data. We then renormalize each year’s data to the 
nominal amplitude and iterate our double-exponential fit. 
This produces a small refinement in our fit. The results 
are shown in Fig. 5. 


In the example of the Daytime Sextantids, the data and 
our derived activity profile differ significantly from our 
previously used profile. In the next section, we summa- 
rize our data and implications for the annual forecast. 
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Figure 4. A double exponential activity profile (black 
line) fit to detrended fluxes (points) corresponding to the 
DSX radiant as measured by CMOR. After fitting is com- 
plete, we remove points that lie more than 50 away from 
the fitted profile — in this case, five such outliers have been 
removed. 


3.3. Revised list of showers 


We were able to characterize the activity profiles for 12 
major meteor showers using CMOR flux measurements. 
In some cases, such as the Orionids (ORI), our fitted pa- 
rameters closely resembled those of [8]. In other cases, 
such as the Daytime Arietids, our activity profile param- 
eters changed substantially. Table 2 shows the list of all 
updated showers, including previous and updated values 
of all parameters. 


The updated showers include the Geminids (GEM); we 
altered the shape of this shower to be much broader. We 
made this choice because CMOR measures fluxes at par- 
ticle sizes that are close to the limiting particles sizes for 
typical spacecraft risk assessments. However, the shape 
of the Geminid activity profile is believed vary with lim- 
iting particle size; we plan to incorporate this behavior 
into the forecast at a future date. 


In many cases, we were unable to improve the activity 
profiles (see Table 3). Showers fell into this category for 
a variety of reasons: their activity may have been too low 
to be clearly detected by CMOR, they may have been 
located in the southern sky and out of CMOR’s range 
(SSG), they may have had brief, inconsistent outbursts 
(DRA), or they may have had such an extended activity 
profile that we risked removing shower flux along with 
the background (NTA, STA, BTA). We hope to improve 
our characterization of many of these showers in the fu- 
ture using flux measurements from other sources. 


Finally, we removed a number of low-activity meteor 
showers from our list. We removed showers if they 
had either few detections in the literature and/or non- 
detections in recent surveys or were removed from the 
IAU meteor shower list. The removed showers (and their 
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Figure 5. A double exponential activity profile (black 
line) fit to detrended fluxes (points) corresponding to the 
DSX radiant as measured by CMOR. The amplitude has 
been permitted to vary by year; in this plot, each year’s 
data has been renormalized to match the average flux 
amplitude. The activity profile used in previous meteor 
shower forecasts (thick gray line) is included for the sake 
of comparison. 


common abbreviations) include the w Scorpiids (OSC), 
e Geminids (EGE or Eps), Comae Berenicids (COM or 
CBe),  Virginids (DLI or mVi), a Centaurids (ACE), 
o Centaurids (oCe), y Velids (GVE), 6 Cancrids (DCA), 
a Carinids (aCa), 6 Leonids (DLE), 6 Pavonids (DPA), 
tT Cetids (TCT or Cet), 7 Aquariids (TAQ), o Cygnids 
(OCY), Southern « Aquariids (SIA or iAZ), Northern 7 
Aquariids (NIA or iAN), Northern 6 Aquariids (NDA or 
dAN), y Doradids (GDO), « Aquariids (KAQ), Piscids 
(PSC or Pis), ¢ Puppids (ZPU), a Monocerotids (AMO), 
x Orionids (XOR or cOr), and Phoenicids (PHO). 


Tables 2 and 3 provide the parameters for our standard 
shower list; this list reflects typical activity for each 
shower. Each year, we adjust those showers that are pre- 
dicted to show unusual activity. For instance, the 2014 
shower forecast included the May Camelopardalids, a 
new meteor shower that was predicted! to make its first 
appearance that year [21, 22, 23]. In 2016, the Perseids 
were predicted to reach approximately twice their usual 
rates and to exhibit multiple peaks in activity”; again, 
these predictions were incorporated into the 2016 annual 
shower forecast. 


4. APPLICATION TO RISK ASSESSMENTS 


The outputs of the meteor shower forecast are designed 
to support quick risk assessments for a spacecraft in low 
Earth orbit: all calculations are done at an altitude of 400 


"http://www.imo.net/potentially-big-meteor-outburst-of-209p- 
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km above the Earth’s surface. We provide fluxes, flu- 
ences, and flux enhancement factors for the worst-case 
scenario in which a spacecraft surface directly faces a me- 
teor shower radiant with no planetary shielding. Because 
hypervelocity impact damage is more closely related to 
kinetic energy than mass, our values are all kinetic- 
energy-limited. Table 1 provides the particles sizes and 
masses that correspond to our limiting kinetic energies 
for a particle velocity of 20 km s~?. 


Figs. 6-9 provide sample outputs from a meteor shower 
forecast that corresponds to typical annual meteor shower 
activity. Fig. 6 shows the zenithal hourly rate over the 
course of the year; while ZHR is not a useful quantity for 
risk assessments, it does illustrate how visual meteor rates 
vary with time. Fig. 7 presents the shower flux, which is 
more applicable to risk assessments. Fig. 7 also includes 
a horizontal line marking the time-averaged total mete- 
oroid flux at 400 km altitude. Note that at smaller par- 
ticle sizes, both the shower flux and total flux are larger, 
but showers contribute a smaller fraction. At the largest 
sizes, the flux is much lower, but showers are relatively 
more significant. The particle size at which showers con- 
stitute half of the total flux is larger than those modeled 
here: 10.5 g at 20 kms7!, or 2.1 MJ (503 g TNT equiv- 
alent). For a density of 1000 kg m~°, this corresponds to 
a particle diameter of 2.7 cm. 


Fig. 8 presents the meteoroid fluence as a function of 
time. In this case, the fluence corresponds to the inte- 
gral of the flux over a 7-hour period. We compute 7-hour 
fluences for our annual forecast, but provide custom fore- 
casts in which the integration period can be specified by 
the customer. Finally, Fig. 9 shows the flux enhancement 
produced by showers over the course of the year. We re- 
port the flux enhancement on a surface facing the meteor 
shower; the contribution of meteor showers to the envi- 
ronment in general is one-fourth as large. 


Our forecast enables quick assessments of the additional 
impact risk posed by a meteor shower. It does not con- 
tain detailed speed and directionality information like that 
provided by MEM [1, 2]. However, we do provide ra- 
diants for showers with high fluxes as part of our an- 
nual forecast. If a quick assessment using the forecasted 
shower flux enhancement indicates that the risk may be 
close to the acceptable threshold, further analysis should 
be done. This analysis should include the calculation of 
the aberrated radiant (taking gravitational focusing and 
the spacecraft’s velocity vector into account), the total 
area presented to this aberrated radiant, and the angle(s) 
of impact onto that surface. The gravitational focusing 
and shielding factors we have used here are averaged over 
a sphere with a radius of Re + 400 km; however, grav- 
itational focusing and shielding can exceed this average 
at some points on this sphere [7], and this, too, should 
be taken into account if the initial risk assessment is con- 
cerning. 


We make one final, cautionary note regarding the use of 
these shower fluxes. We have endeavored to provide a 
worst-case scenario for the increase in particle flux. How- 


ever, the width and depth of an impact crater can also de- 
pend on the impact angle. Thus, if a spacecraft surface 
faces the shower radiant, and the ballistic limit equation 
relevant to that surface predicts greater damage when the 
impact is normal (as in the modified Cour-Palais BLE for 
aluminum, [18]), the risk could be elevated by an even 
greater factor than the forecast is able to describe. 
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Table 3. Forecasting parameters for 17 meteor showers for which our activity profile parameters (Xo, Bp, and Bm) 
remained unchanged. However, we have in some cases updated the shower name and corresponding three-letter code to 
comport with the International Astronomical Union’s meteor shower list. Finally, we include peak zenithal hourly rate 
(ZHRg), the shower’s velocity at the top of the atmosphere (V0 km), and the population index (r). 
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Figure 6. MEO annual meteor shower forecast sample output. This plot shows the total zenithal hourly rate (ZHR) due to 


meteor showers over the course of one year. Peaks corresponding to major meteor showers are labeled with the shower 
code. 
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Figure 7. MEO annual meteor shower forecast sample output. This plot shows the total particle flux due to meteor showers 
over the course of one year for four particle energies (see Table 1). Peaks corresponding to major meteor showers are 
labeled with the shower code. The total flux (sporadic and shower combined) corresponding to each particle energy is 
also shown as a horizontal line. Note that for small particle sizes, meteor showers do not exceed the average flux, while 
at large sizes, many showers exceed the average. A similar plot is presented in Fig. 5.2 of [7]. 
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Figure 8. MEO annual meteor shower forecast sample output. This plot shows the 7-hour particle fluence (i.e., the flux 
integrated over a 7-hour period) due to meteor showers over the course of one year for the smallest reference particle 
size. Peaks corresponding to major meteor showers are labeled with the shower code. 
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Figure 9. MEO annual meteor shower forecast sample output. This plot shows the enhancement of the total flux due to 
meteor showers for the smallest two reference particle sizes. Peaks corresponding to major meteor showers are labeled 
with the shower code. Fluxes, fluences, and enhancement factors assume that the spacecraft surface in question is oriented 
to face the shower radiant, maximizing the intercepted flux. 


